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Truncated Levy flights are stochastic processes which display a crossover from a heavy-tailed Levy 
behavior to a faster decaying probability distribution function (pdf). Putting less weight on long 
flights overcomes the divergence of the Levy distribution second moment. We introduce a fractional 
generalization of the diffusion equation, whose solution defines a process in which a Levy flight of 
exponent a is truncated by a power-law of exponent 5 — a. A closed form for the characteristic 
function of the process is derived. The pdf of the displacement slowly converges to a Gaussian in 
its central part showing however a power law far tail. Possible applications are discussed. 

PACS numbers: 02.50.-r; 05.40.Fb 
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Since Levy flights have been introduced into statistical 
physics, it has become clear that special attention must 
be given to the fact that due to their heavy tails they are 
characterized by diverging moments. A few approaches 
have been suggested to overcome this divergence. These 
include the introduction of the concept of Levywalks [J , 
confining Levy flights by external potentials and in- 
troduction of truncation procedures 0, . Each of the 
approaches represents a different physical situation, but 
they all made it possible for Levy processes to be ap- 
plicable in a variety of areas ranging from Hamiltonian 
dynamics i5| to spectral diffusion in single molecule spec- 
troscopy Iq], from bacterial motion 7] to the albatross 
flights <8| and also in analysis of economical data |^ ^3 ■ 
Here we concentrate on the truncation of the flights. 

In many cases the Levy-flight behavior corresponds to 
intermediate asymptotics. At very large values of the 
variable some cutoff enters, so that the moments exist. 
Truncated Levy flights, a process showing a slow con- 
vergence to a Gaussian, were introduced by Mantegna 
and Stanley and have been since used especially in 
econophysics, see Refs. d [13. The truncated Levy 
flight is a Markovian jump process, with the length of 
jumps showing a power-law behavior up to some large 
scale. At larger scales the power-law behavior crosses 
over to a faster decay, so that the second moment of 
the jump lengths exists. In this case the central limit 
theorem applies, so that at very long times the distri- 
bution of displacements converges to a Gaussian; this 
convergence however may be extremely slow. The orig- 
inal work concentrated on numerical simulations of the 
process which assumed a 0-function cutoff. Koponen Q 
slightly changed the model by replacing the (^-function 
cutoff by an exponential one and obtained a useful an- 
alytical representation for the model. However, further 
investigations have shown that the models with sharp [9- 
function or exponential) cutoffs predicting a Gaussian or 



an exponential tail of the pdf are not always appropriate 

In mathematical physics, it is often convenient to have 
a deterministic equation for the pdf of a process, an ana- 
logue of the diffusion or Fokker-Planck equation (FPE), 
to be solved under given initial and boundary conditions. 
For the case of anomalous transport, fractional general- 
izations of such equations may be relevant |ll| . How- 
ever, these classes of equations are valid only for processes 
showing exact scaling in the force-free limit. The trun- 
cated Levy flights are not such a process. As mentioned, 
in the course of time the process displays a crossover 
between the two regimes: at shorter times the character- 
istics of this random process behave as those of a Levy 
flight, while at long times they are close to ones for nor- 
mal diffusion. Processes showing a crossover can often be 
described by equations containing derivatives of different 
order in the same variable. A known example here is 
the telegrapher's equation with the second and the first- 
order temporal derivatives, describing the crossover from 
ballistic transport to the diffusion behavior. However, 
fractional generalizations of the telegrapher's equation 
describe Levy-walk-like processes 12] and not truncated 
Levy flights. 

The equation we propose for truncated Levy flights 
with the power-law cutoff has the following form: 



dp{x,t) _ d^p{x,t) 



d\x\ 



dt 



= D 



(1) 



where D is the diffusion coefficient governing the long- 
time asymptotic behavior, and the scale factor Ca = 
D/Ka is a coefficient governing the intermediate-time 
Levy-like one. The dimension of Ca is [Ca] = \L'^~°'\- In 
Eq-O g^iQ denotes the symmetric Riesz-Weyl operator 
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|l3| . which can be expressed through 



— [^/W = 



d\xf" - ' 2cos(^/3/2) 
for < /? < 2, /3 7^ 1 and 



dlxf 



/(^) 



da; 



(2) 



(3) 



for /3 = 1, where -oaD^ and xD_^ are the corre- 
sponding Riemann-Liouville operators, and H denotes 
the Hilbert transform 



H(f>=- 



(4) 



The operator defined by Eqs. l(2Jl and Q is a 
fractional generahzation of the second derivative: in 



Fourier- representation 



(x) simply corresponds to 



-\k\ 4>{k), which for /3 = 2 gives us a known form 
-k'^(j>{k). This is the reason to include the "minus" sign 

(a;) corresponds to 



9^ 



in Eq.lO. Note however that 

Eq.lQJ is a special case of the distributed-order diffu- 
sion equation: 
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o\x\ 



dp{x, t) 



dt 



corresponding to a weight function /(a') = S{2 — a') + 
6 {a — a') which describes a crossover between a — 2 
(normal diffusion) and to a < 2 (Levy-like superdiffu- 
sion). It can be shown, that contrary to the equations 
used in 0,^3 which describe processes getting more and 
more anomalous in the course of time (retarding subd- 
iffusion or accelerating superdiffusion) , equation de- 
scribes processes getting less anomalous, such as in our 
case, tending to normal diffusion. General properties of 
this equation will be considered elsewhere. 

Using the Fourier-representation of the Riesz-Weyl 
derivative in Eq. ^ we get for the characteristic function 
of the distribution: 



l + Co.\k\' 



dfjkA) 
dt 



= -Def{k,t). 



(6) 



The Green's function of this equation, corresponding to 
the initial condition /(fc,0) = 1 (i.e. p{x,0) — S{x)) then 
reads: 



The function f{k, t) is differentiable twice for each a; 
its second derivative /"(fc,i)|fc^o ~ —2Dt, so that the 
second moment of the distribution evolves in a diffusive 
manner: ^x^(<)) = 2Dt, as in normal diffusion. How- 
ever, in the intermediate domain of x the distribution 
shows the behavior typical for Levy flights; namely for k 
large enough, i.e. for Cq. ^ 1, the characteristic 

function has the form 



/(fc,t) = exp (^-^|fc|"t^, 



(8) 



i.e. corresponds to the characteristic function of the Levy 
distribution. Assuming a < 2 we get the following ex- 
pansion for f{k, t) near fc = 0: 



f(k,t) ~ I- Dtk^ + DCo,t\k\ 



(9) 



From this expression it is evident that /(fc, t) always lacks 
the fourth derivative at fc = (for 1 < a < 2 it even lacks 
the third derivative), which means that the fourth mo- 
ment of the corresponding distribution diverges. The ab- 
sence of higher moments of the distribution explains the 
particular nature of the truncation implied by our model: 
The Levy distribution is truncated by a power law with a 
power between 3 and 5. Thus, for all < a < 2 the corre- 
sponding distributions have a finite second moment and, 
according to the central limit theorem (slowly!) converge 
to a Gaussian. For the case < a < 1 (when ( | ) < oo) 
the speed of this convergence is given by the Berry-Esseen 
theorem, as noted in Ref.jl6l|. The convergence criteria 
for 1 < a < 2 can be obtained using theorems of Ch. 
XVIofRef. [13. 

This transition from the initial Levy-like distribution 
to a Gaussian is illustrated in Fig.l, obtained by a numer- 
ical inverse Fourier-transform of the characteristic func- 
tion, Eq.©. Here the case a = l,Z3 = C= lis shown. 
To put the functions for t = 0.001 and for t = 1000 on 
the same plot we rescale them in such a way that the 
characteristic width of the distribution W{t) (defined by 

/o^*^*^ p(a;, t)dx — j) is the same. The behavior of the pdf 
to be at the origin p(0, t) as a function of t is shown on 
the double logarithmic scales in Fig. 2. Note the crossover 
from the initially fast decay p(0,t) cx t^^/°' (Levy su- 
perdiffusion) to the final form p(0,i) oc i^^/^ typical for 
diffusion. 

The asymptotics of the pdf at large x is determined by 
the first non-analytical term in the expansion, Eq.Q, i.e. 
hy DC at \ k\ ~". By making the inverse Fourier transfor- 
mation of this term and using the Abel method of sum- 
mation of improper integral, we get 



/(fc,t) = exp 



l + Ca\k\ 



(7) 



p{x,t) 



r(5-a) sin(7ra/2) DC^t 



v,5 — a 



(10) 



We postpone the proof of the fact that this is indeed a 
characteristic function of some probability distribution 
p(a;, t) until later on and discuss first the main properties 
of such a solution. 



where we use that Jp°° cos^ = sin(7rQ;/2)r(5 — a). 

Thus, in our case the Levy-distribution is truncated not 
by a 9- or exponential function, but by a steeper power- 
law, with a power = 5 — a. 
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FIG. 1: The rescaled pdf p(^) — W{t)p{x,t) is shown for 
t = 0.001 (dotted hme) and for t — 1000 (dashed line) as a 
function of a rescaled displacement ^ — x/W{t). The corre- 
sponding thin lines denote the limiting Cauchy and Gaussian 
distributions under the same rescaling. 



For example a Levy flight truncated by another, faster 
decaying power-law, is a much better model for the be- 
havior of commodity prices. Thus, the discussion in 
Ref.l^ shows that the cumulative distribution function 
of cotton prices may correspond to a power-law behavior 
of 1 — F{x) = p{x)dx oc a;^" with the power a — 1.7 
in its middle part and with the far tail decaying as a 
power-law 1 — F{x) cx x~^ with /3 « 3. Thus, our equa- 
tion (which is definitely the simplest form of the equation 
for truncated Levy flights) adequately describes this very 
interesting case giving f3 ~ 3.3. It is highly probable that 
fractional equations of the type considered here might be 
a valuable tool in economic research. 

Let us now prove that the solution p{x, t) is a pdf, i.e. 
a non-negative normalized function of x for any t. The 
normalization is trivial and follows from the fact that 
/(O, t) — 1 for all t. Let us now prove the non-negativity 
of the solution. 

We start from defining a function G{u,t), u > 0, such 
that its Laplace transform in variable u is 

G{s,t) = J^^du e-^-Giu,t) = exp ^-—^f-^t^ 

(11) 

with Aa = Ca/D^^"/^. By comparing Eqs.UIll and 10, 
one sees that the characteristic function f{k,t) can be 
rewritten in the following form: 
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FIG. 2; Shown is the pdf to be at the origin p(0, t) as a func- 
tion of time, see text for details. Note the double logarithmic 
scales. The dashed line has the slope —1, and corresponds 
to the superdiffusive decay; the dotted line has the slope of 
— 1/2, as in the case of normal diffusion. 



j{k,t)= e-''°'''G{u,t)du. (12) 
Jo 

Indeed, the transition from Eq.Q to Eas. ll2|l is noth- 
ing else but the change of variable s — > Dk^ in Ea. (|ll|l . 
It is clear that G(0, t) = 1 and, moreover, as we pro- 
ceed to show, G{s, t) is completely monotonic, i.e. it is 
non-negative, and the signs of its derivatives alternate. 
Then, according to Bernstein's theorem G'(s,i) is a 
Laplace-transform of some probability density. Now, we 
can perform the inverse Fourier-transform in the Eq. I|12|) 
and get 

Jo VAirDu \ 4:DuJ 

which is a nonnegative function (since the integrand 
is a product of two non-negative functions). Ea. H13|l 
provides a subordination transformation: the truncated 
Levy flights can be considered as a process subordinated 
to a Wiener process under the operational time given by 
the function G{u,t) The small-u behavior of this 

function (corresponding to the large-s one of its Laplace- 
transform, G{s,t) = exp (— A~^s"''^i) ) is approximately 
a one-sided (extreme) Levy law of index a/2. However, 
at large u this Levy law is truncated. 

We now give a proof that the function G{s, t) — 

exp ^— is a completely monotonic function 

in variable s. This function has a form exp(— 7/;(s)) and 
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therefore is completely monotonia if the function ip{s) 
is positive and possesses a completely monotonic deriva- 
tive. In our case ip{s) — s/(l + As^) with b = 1 — a/2, 
< 6 < 1, A > 0. Its derivative is given by 



1 



1-6 



As" 



1 + Asb 



(14) 



This function is a product of two functions. The first one 
is completely monotonic since it has a form g{h{s)) with 
g{y) = 1/(1 + y) being completely monotonic and with 
h{s) — As" being a positive function with a completely 
monotonic derivative. The second function has the same 
form, now with g{y) = 1 — by/{l + y). This function 
g{y) is positive for all ?; > 0, and its derivatives read 
c,(")(y) = (-l)»6n!(l + yy-\ 

The subordination property also sheds light on the pos- 
sible nature of truncated Levy distributions in economic 
processes. The truncated Levy process can be interpreted 
as a simple random walk with a finite variance. However, 
the number of steps of the random walk (the number of 
transactions) per unit time is not fixed, but fluctuates 



strongly. The implications to economics of such mod- 
els were considered in [l9l |. In our case the distribution 
function of the number of steps has itself a form of a 
truncated one-sided Levy law. 

Let us summarize our findings. We proposed a frac- 
tional generalization of a diffusion equation which de- 
scribes power-law truncated Levy flights, a random pro- 
cess showing a slow convergence to a Gaussian. We show 
that the solution of this equation is a pdf and give nu- 
merical results for the case a — 1. Moreover, we argue 
that the truncated Levy flights can be represented as a 
random process subordinated to a Wiener process, which 
might be helpful in econophysical applications. We end 
by noting that the equation discussed is a special case of 
distributed-order fractional diffusion equations. Modifi- 
cations of our equation should be able to describe other 
types of truncation; however, the equation discussed here 
is definitely the simplest one. 
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